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ABSTRACT 

Pulsar timing arrays (PTAs) are being used to search for very low frequency gravitational 
waves. A gravitational wave signal appears in pulsar timing residuals through two compo- 
nents: one independent of and one dependent on the pulsar's distance, called the 'Earth term' 
(ET) and 'pulsar term' (PT), respectively. The signal of a burst (or transient) gravitational 
wave source in pulsars' residuals will in general have the Earth and pulsar terms separated by 
times of the order of the time of flight from the pulsar to the Earth. Therefore, both terms are 
not observable over a realistic observation span, but the ETs observed in many pulsars should 
be correlated. We show that pairs (or more) of pulsars can be aligned in such a way that the 
PTs caused by a source at certain sky locations can arrive at Earth within a time window short 
enough to be captured during a realistic observation span. We find that for the pulsars within 
the International Pulsar Timing Array (IPTA) ~ 67 per cent of the sky produces such align- 
ments for pulsars terms separated by less than 10 years. We compare estimates of the source 
event rate that would be required to observe one signal in the IPTA if searching for the cor- 
related ETs, or in searching via the PTs, and find that event rates would need to be about two 
orders of magnitude higher to observe an event with the PTs than the ETs. We also find that 
an array of hundreds of thousands of pulsars would be required to achieve similar numbers of 
observable events in PT or ET searches. This disfavours PTs being used for all-sky searches, 
but they could potentially be used target specific sources and be complementary to ET only 
searches. 
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1 INTRODUCTION 

In the late 1970s it was first suggested that precision timing 
of pulsars co uld be used to detect very low frequency gravita- 
tional waves jSazhin|[l978l : |Petweilej[l979l) . This lead to early 
searches for a cosmol ogical stochastic backgrou n d of nanoHz grav- 
itational waves (e.g. iHellings & Downs! 1983 ; Romani & Tavloj 
1983; lDavise~a l. 1985; R awlev et al.ll 19871 : IStinebring et al J 199(3 : 
Kaspi. Taylor & Rvballl994l) . The first attempt to construct an ar- 



ray of pulsars for use as a gr avitational wave detecto r (amongst 
other applications) was that of iFoster & Backej dl99Qh . More re- 
cently sev eral groups around the world (the Parkes Pulsar Tim- 
ing Array ( Verbiest et al .1120 lOl) , the European Pulsar Timing Array 
dFerdman et alj|2010h and the N orth American Nanohertz Obser- 
vatory for Gravitational Waves dJenet et al J 120091 : iDemorest eUd] 
l2012h ) have worked to set up and perform precision timing of a se- 
lection of stable millisecond pulsars. The aim of these is to detect 
low frequency gravitational waves from objects such as supermas- 
sive binary black holes (SBBH). These are now being combined 
into a concerted worl d-wide effort to for m an International Pulsar 
Timing Array (IPTA) faobbs et alEoioh . 
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Initial gravitational wave searches using pulsar timing focused 
on looking for a cosm ological stochastic ba ckgr ound. Following 
the the oretical work of Ijaffe & Backed d2003h and lWvithe & Loebl 
d2003h the focus has more recently shifted to findi ng a s t ochas- 
tic background from multiple SBBHs. iLommen & Backed J200ll) 
performed the first searches for individual quasi-monochromatic 
SBB H sources (with the more recent theoretical work of 
Sesana, Vecchio & Volonteri 2009, providing more impetus for 
this), with this method being used to rule out a putative electro- 
magnetic observa tion of such a system in the radio galaxy 3C 66B 
dJenet et al.ll2004r) . Now th ere are many proposed methods to de- 
tect stochastic sources fe.g.lMcHugh et alJI 19961 : Ijenet et al.|[2005l : 
lAnholm et aDl2009l : IvanHaasteren et alJl2009h . and several recent 
searches have been performed providing limits on their emis- 
sion (e.g.lJenet et al.l2006l:lYardlev et al.l201ll ; Ivan Haasteren et all 
l201ll : IDemorest et alj |2012). There are now also many proposed 
methods to search f or individual quasi-mono c hromatic (or contin- 
uous) sources ( e.g. ICorbin & Cornish) l20ld: lYardlev et al.ll2010l: 
Lee etalJ koil: Babak & Sesana 200: lEUis. Jenet & McLaughlinl 



duration transient (or burst) sources (e.g. 


van Haasteren & Levinl 


2010tlPshirkov Baskaran & PostnovfeOlCt 


Finn & Lommenll2010l: 


Cordes & Jenet 2012|). A review of the many gravitational wave 
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search avenues currently being explored can be found in lLommenl 
d2012h . 

Examples of burst sources could be the final in spiral or 
parabolic encoun ters of SBBHs jFinn & Lommen l201Qh. or cos- 
mic s t ring cusps jLeblond. Shlaer & Siemensll2009 ; lBinetruv et alj 
120091 : iKev & Cornish! 120091) . Unlike the bur sts searched for i n 
ground-based gravitational wave searches (e.g. I Abbott et al 1l2QQ9h . 
which are generally classed as events lasting of the order of 
milliseconds-seconds, these bursts would last from months to 
years, but importantly they are transient rather than continuous sig- 
nals. 



1.1 Signals in pulsar timing arrays 

Pulsar observers measure the time of arrival (TOA) of pulses, which 
can be thought of as ticks of a clock. Timing residuals are the differ- 
ence between the observed TOAs and a best fit model of the time of 
arrival that is dependent on many parameters such as the pulsar fre- 
quency, frequency derivatives, sky position and any binary system 
parameters if appropriate. Any unmodelled components, such as a 
potential gravitational wave signal, would remain in the residuals. 

A gravitational wave signal appears in pulsar timing residuals 
through two terms: an 'Earth term' and a 'pulsar term', which we 
shall refer to as ET and PT from now on. The ET is independent 
of the pulsar distance and the PT is delayed in time from the ET 
by an amount proportional to the pulsar distance. This 'two-pulse' 
response is more generally the case for an y signal in a single-arm 
one-w ay gr avitational wave d etectoj]] e.g. lEstabrook & Wahlquistl 
dl975h and lDetweileil dl979h . In this case the pulsar and the re- 
ceiver on the Earth represent the ends of the arms of the detector. A 
signal in the ET will be simultaneous in all observed pulsar residu- 
als, whereas the PT will be delayed by 



At = (1 + k • h)d/c, 



(1) 



where k is a unit vector along the gravitational wave propaga- 
tion direction (pointing from the source to the Earth), n is a unit 
vector pointing along the Earth-pulsar line-of- sight, and d is the 
distance to the pulsar. If we define a plane perpendicular to k 
such that the Earth is at the origin and k — {0,0,-1}, and 
n = {sin cos 0, sin sin 0, cos 0}, where is the polar angle 
between the source and pulsar, and (j) is the azimuthal angle from 
some arbitrary point (which for the moment does not matter), then 

At = (l-cos0)d/c. (2) 

As stated in lFinn & Lommenl J2010L unless the pulsar and source 
are closely aligned on the sky, this time delay will be large. So, 
for a transient burst source if the ET is observed then the PT 
will not appear in the residuals over the typical period of pul- 
sar observations of one-to-two decades. For this reason the PT 
is often ignored. The inclusion of the PT, and extra information 
that can be gathered from it (e.g. pulsar distance measurements 



dSazhinll97 8l;lLommenl200ll;|jenet et al . 2004; Yar dlev et al.120101 : 
>Kj;" 



ICorbin & Cornishll2010l: iLee et alJbOllh T source distance via par- 
allax foeng & F inn 2011 ) or studyi ng a signal at different times in 
its evolution jMingarelli et al.ll2012h \ has been studied with regard 
to continuous gravitational wave signals, but not for burst sources. 



1 Note that for ground-based gravitational wave detectors, where the wave- 
length of the gravitational wave is far longer than the arm length, the 'two- 
pulse' struc ture and fre quency dependence disappear from the detector re- 
sponse (e.g. En32009b . 



However, pulsar timing arrays (PTAs) consist of many pul- 
sars. This brings the possibility that among pairs, or more, of pul- 
sars there might be fortuitous lines-of- sight where the time delay 
between PTs in particular pulsars is small enough to be within an 
observational data set. This is the possibility we explore in Sec- 
tion |3 of this paper. 

In a standard burst search using just the ETs all pulsars in the 
array can be coherently analysed giving the maximum possible sen- 
sitivity, but temporally it is limited to signals that are passing the 
Earth now (or at least to within the decade or so over which the 
PTA is observing). Looking for signals in the PT will not give as 
good sensitivity as that from a coherent ET search, because only 
two (or a few) pulsars can be coherently combined, but each pul- 
sar pair will provide a different time baseline in which to search. 
These signals would otherwise not be observed at all in an ET-on\y 
search. So, given the temporal coverage it provides it could still be 
a worthwhile area to search. We will explore this more in Section[3] 



2 INTER-PULSAR TIME DELAYS 

In this section we will con sider the 30 pulsars of the IPTA as given 
in lFinn & Lommenl j2010h . Here we look at the minimum time de- 
lay between the PTs for all pairs of pulsars in the array for sources 
located across the sky. For a pulsar i with right ascension and decli- 
nation on and Si, distance d\, and a source position a s and S s , with 
a distance L emitting at to, the time of arrival of the PT at the Earth 
will be 



U = to + L/c + (1 - cos 6i)di/c, 
where the angular separation between the source and pulsar is 



(3) 



0i = arccos {cos (a) cos (b) + sin (a) sin (b) cos (\cti — ot s \)}, 

(4) 

with a = 7r/2 — Si and b — tt/2 — S s . So, for a pair of pulsars i and 
j the difference between the time of arrival of the PT for a source 
at a given sky position is 



AU 



\U-tj \ = |(1 - cos 0i)di - (1 - cos 0j)dj |/c. (5) 



A schematic of this set-up for a pair of pulsars is shown in Fig- 
ure Q] It is important to note that the delay St in the figure is not the 
time delay between the PTs, because the information that the grav- 
itational wave has influenced the pulsar still has to travel to Earth 
encoded in the electromagnetic pulses. Thus the real time differ- 
ence is given by Afoi- Figure |2] shows how the time delay, Afoi, 
between the times-of-arrival of the PTs observed in pulsar residuals 
at Earth for a PTA consisting of a pair of pulsars (J1455-3330 with 
a = 14 h 55 m 47!9, S = -33 o 30'46"3 and distance of 0.74 kpc, 
and J2129-5721 with a = 21 h 29 m 22"7, S = -57°21 / 14 / /l and a 
distance of 0.53 kpc) changes for sources located across the whole 
sky. It can be seen that the minimum in the delay forms a ring on the 
sky. The overlap of many such rings can be seen later in Figure |4] 

If we choose a constraint such that in Equation \5\ AUj < 
10 years then Figure [3] shows the valid areas of the 0i, 0j, di and 
dj parameter space that fulfil this criterion. It can be seen that in 
the extreme cases where one pulsar is closely aligned on the sky 
with the source, and the other is on the opposite side of the sky to 
the source (e.g. Oij « 0° and 0j,i ~ 180°) then the pulsar distance 
ratio must be very large to fulfil the criterion. 

In this paper we will assume a time delay between PTs of 
less than 10 years is reasonable for searches on current, or near 
future, sets of data. Using a grid of source sky positions and all 
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Figure 1. A schematic diagram showing the geometry of two pulsars, a 
passing plane gravitational wave and the Earth 0. From top to bottom the 
planes correspond to the gravitational wave intersecting the first pulsar, then 
the second pulsar and then the Earth. The time delay between the PTs is 
given as A £21- 




Figure 2. The time delay between the PTs from Equation \5\ as a func- 
tion of source sky position for a pulsar pair consisting of J1455-3330 
(a = 14 h 55 m 47!9, 5 = -33°30 / 46 / /3 and distance of 0.74 kpc) and 
J2129-5721 (a = 21 h 29 m 22!7, S = -57°21 / 14 / /1 and a distance of 
0.53 kpc). 



the pulsars in the IPTA, we have calculated the minimum value 
of At 21 in years for all pulsar pairs, the sky area for which this 
delay is less than 10 years and the total number of pulsar pairs for 
which this delay is less than 10 years. These are shown as Hammer 
projections onto the sky in Figured Obviously for longer spans 
of data longer time delays can be contemplated and more pulsar 
pairs will be usable. Figure [4(b)] shows the sky area available if we 
exclude all sky positions for which no Prtime delay is less than 10 
years. This shows that 67 per cent of the sky will contain at least 
one pair of pulsars with a time delay between PTs of less than 10 
years. Residuals for any such pair could then be cross correlated 
with the appropriate delay applied for the sky position. 

Figure [4(c)] shows the total number of pulsar pairs that would 
have time delays of less than 10 years for each sky position of the 



/ H ci 
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Figure 3. The sheets represent regions in the 



, di and dj parame- 



ter space for a pair of pulsars that fulfil the criterion from Equation \5\ that 
Atij < 10 years. 



gravitational wave source. We see that there are some small patches 
of the sky where up to 11 pulsar pairs are usable. However, it is 
unlikely that multiple pairs of pulsars would have the same indi- 
vidual signal in them for a given sky position. Unfortunately this 
means that generally coherent analyses between multiple pairs are 
not viable, but having multiple pairs gives you a higher effective 
temporal observation span (as discussed in Section [3]). There are 
sky locations for which more than two IPTA pulsars are aligned 
as such to fulfil the time delay criterion, but the sky area covered 
is general small. Figure [5(a)! shows that for the IPTA we find 5.6 
per cent of the sky for which the PT for three pulsars are within 
10 years of each other, 0.3 per cent for which there are four pulsars 
and 0.005 per cent for which there are the maximum of five pulsars. 



2.1 Expanding the PTA 

In the future when more pulsars are added to an array more sky 
areas with multiple pulsar overlap (i.e. > 2) may be available, 
in particular with the addition of many pulsars within the same 
globular cluster, which are physically separated by short (of or- 
der years) delays. An example of how co-located pulsars in glob- 
ula r clusters can be used in gra v itation al wave detection is given 
by Ijenet Creighton & Lommenl d2005h . As a first step at testing 
how increased numbers of pulsars improve prospects, we have cre- 
ated a fake array by taking all millisecond pulsars currently in the 
Australian Telescope N ational Facility (ATNF) pulsar catalogue 
(Manchester et al.ll2005l) with frequencies greater than 50 Hz, but 
not associated with a globular cluster, leaving 90 pulsars. These are 
just used as an example of what happens when adding more pul- 
sars, however it is not known, or expected, that it will be possible 
for some, or all, of these specific pulsars to be included in a future 
array. With this array, and assuming a well-known distance to these 
pulsars, the source locations for which at least a pair of pulsars has 
PTs within 10 years show over 99.8 per cent sky coverage. We find 
that the number of pulsars for which the PTs are separated by less 
than 10 years can be up to 8, although this is for one unique sky 
position (a = 18 h 32 m 23!0 and 6 = 01°34 / 44 / /21) (given a sky 
pixel of area ~ 0.7 square degrees). The overall number of pulsars 
with PTs within 10 years for all source locations across the sky can 
be seen in Figure [5(b)] In terms of sky area there is 42 per cent of 
the sky for which 3 pulsars have PTs within 10 years, 7 per cent 
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(c) 

Figure 4. |(a)| A Hammer projection sky map giving the minimum time delay between PTs over all IPTA pulsar pairs for a source located at each position in 
the sky. The colour range gives the time delay in years, [(b)] A map showing only those parts of the sky (black) where a source will produce a signal separated 
by less than 10 years within at least one pulsar pair, [(c)] A map showing the number of pulsar pairs for which a source at that sky position would produce a PT 
time delay of less than 10 years. All maps show crosses for the locations of the IPTA pulsars. 
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with 4 pulsars, ~ 1 per cent with 5 pulsars, ~ 0.08 per cent with 
6 pulsars and ^0.01 per cent with 7 pulsars. More pulsars in the 
array can therefore increase the sensitivity for certain sky areas. It 
should be noted that the sky coverage may be biased by selection 
effects of observable stable, well-timed, pulsars. 



2.2 Distance uncertainties 

The initial discussions assume that we know the distances to all the 
pulsars perfectly, but in reality there will be uncertainties on these 
distances. Currently these are optimistically of order 10 per cent, 
but by the time of the Square Kilometre Array (SKA) errors may 
be reduced to less than 1 per cent for the vast majority of IPTA 
pulsars, with some of the closest known to maybe a few tenths of a 
percent dSmits et alJuOlll) . This uncertainty in the distances means 
there could be some pulsar pairs that actually lie outside of the 
required PT delay criterion. Conversely, accurate distance estimates 
could move some pulsar pairs to within the PT delay criterion. On 
average it would be expected that as many pairs move within the 
criterion as move out of it, but unless the distance uncertainties are 
included in a search it is not known which move in or out. This 
means that in reality the area covered in Figure |4(b)| would look 
different, although the total area covered would be approximately 
the same. 

The effect of this means that for a real search the whole sky 
would have to be covered and the uncertainties in the distances to 
each pulsar taken into account. So, for each sky position every pul- 
sar pair would have to be tested with the distances of each varied 
within the uncertainty range. 



3 SEARCH COMPARISONS 

Many methods have been developed to detect, and estimate param- 
eters, for modelled and unmodelled bursts of gra vitational waves in 
ground-based detector s (e.g. a selection in c lude s Clark et al ]|2007l: 
iKlimenko et all 120081: [Abbott et all [20091 : ISearle. Sutton & Tintd 
2009: ISutton et alJ2O10|). space-based detectors teabak et al.l2010h 
and PTAs (e.g. iFinn & Lommenl l2010l : Ivan Haasteren & Levinl 
[20l3). It would be relatively straightforward to use similar meth- 
ods in a search given the situation that we have presented. So, here 
rather than define another method we will compare aspects of a PT- 
only search with a more standard i^T-only burst search for a simpli- 
fied IPTA. The aspects we compare are their sensitivities, their ef- 
fective observation times, and from these, the estimated event rates 
required to give a detection. 

For this we will compare the signal-to-noise ratio that could 
be recovered for a single burst source detected in the ETs of the 
IPTA with that which could be recovered for a Pr-only search. We 
will make the simplifying assumption that all 30 pulsars within the 
array have residual noise that is white and Gaussian with standard 
deviation of 100 ns, all data spans the same 20 years period and 
they are sampled every 30 days. We will require that to be detected 
a signal must have a signal-to-noise ratio (p) above a threshold of 
10. We will also assume that all the pulsar distances are precisely 
known rather than trying to assess the effects searching over some 
distance uncertainty has on the threshold required for detection^. 



3.1 The signal model 

In these tests we assume a simple signal model of a sine-Gaussian 
burst. This is defined in the timing residuals of a pulsar (indexed by 
i) by 

r t (T) = A+F+(^,fc,ni)[cos(wt? +0 o )exp(- 



- cos (cjtf + 0o ) exp 



2r 2 



)] 



A x F x (^,fc,n,)[sin {ujf % + O ) exp ( - 
-sin« + o )ex P (-^f) 



(6) 



where A+ = A(l + cos\)/2 and A x = Acosl for an ampli- 
tude A and source inclination angle i, r is the Gaussian width, uj is 
the angular frequency, 0o is the phase at the midpoint of the burst, 
and the definitions of the antenna patterns F + / x are described in 
Appendix |A] The times are defined so that the PT is described by 
the terms containing t p — T — U, where T is the pulsar proper 
time at the solar system barycentre and ti is time at the observed 
midpoint of the burst in the PT (as seen for a particular pulsar), 
and therefore the ET time is described by the terms containing 
t e = T — U + (1 + k • fii)di/c. This could be expanded to a 
physical source model such as the parabolic encounter of two su- 
permassive black holes used in IFinn & Lommenl d2010h (which is 
qualitatively similar to a sine-Gaussian), or a more generic burst 
model as employed in ground-based detector searches. 



3.2 Event rates for detection 

It is useful to try and compare the source event rates that would be 
required to detect a signal in the ETs of all pulsars in the IPTA to 
that which would be required to detect a signal in the PTs of pairs of 
pulsars in the IPTA. We have made some general estimates of this 
for the sine-Gaussian signal described in Section 13.11 with fixed 
parameters (following the definition in Equation [6]): A — 500 ns, 
U — 2tt x 4yr _1 , ip = Orads, 0o = 2.0rads, r = 100 days and 
with the burst centred on the midpoint of the observations. 

We define the signal-to-noise ratio to a given signal, r, for a 
set of N pulsars as 



H££- 



1/2 



(7) 



where ra is the number of residual data points, and at is the noise 
standard deviation, for the i th pulsar (working in the time domain). 



3.2.1 Search sensitivity and range 

Using Equation [7] we have calculated the signal-to-noise ratio for 
the sine-Gaussian source defined above at each point in the sky for 
the IPTA pulsars. Figure [6(a)] shows the signal-to-noise ratio that 
this source (if optimally oriented with cos^ = ±1) could be ob- 
served at if located across the sky and seen in the ETs for the entire 
IPTA. The flower-like patterns in the response across the sky come 



2 In reality the signal-to-noise ratio threshold at which to set a detection 
would be lower for the ET search than the PT search. This is because the 
background level (e.g. false alarms due to noise) for a signal being observed 



would be lower due to more detectors being used and the smaller parameter 
space from not having to search over pulsar distances uncertainties. There- 
fore, the horizon distance for ET searches could be increased by a small 
factor of the order of 1.25 reducing required rates by a factor of ~ 2. 
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(b) 

Figure 5. |(a)| The number of pulsars with PTs separated by less that 10 years for given source sky locations using the ITPA. [(b)] The number of pulsars with 
PTs separated by less that 10 years for given source sky locations using a fake PTA of 91 pulsars. 



from the antenna patterns of pairs of pulsars as seen in Figures IaTI 
and lA2l The range of signal-to-noise ratios for this particular source 
over the sky was between 40.7 and 14.8, whereas if the the worst 
case source orientation were chosen (cos l = 0) then the signal-to- 
noise ratio was between 10.0 and 2.5. 

Figure [6(b)] shows the maximum signal-to-noise ratio that this 
source (if optimally oriented with cos i — ±1) could be observed at 
if located across the sky and seen in the PTs of pairs (or multiples 
when applicable) of IPTA pulsars. We define regions of the sky 
with no pulsar pairs PTs separated by less than 10 years as having 
zero signal-to-noise i.e. we assume that we could not see, or would 
ignore as potential noise, a PT signal only observed in one pulsar. 
In this case what we mean by maximum is that for locations for 
which multiple pulsar pairs could contain signals we have taken 
the signal-to-noise ratio of the loudest pair. Over the sky area that 
could be observed (~ 67 per cent of the sky) the range of signal- 
to-noise ratios for the maximum and optimally oriented case was 
between 19.4 and 0.4, whereas for the worst case orientation (and 
with the quietest pulsar pair chosen) the ranges was between 4.6 
and 7x 10~ 4 (see Figure [6(c)]). 



all 30 pulsars) and a PT search (using only a pair of pulsars for each 
sky position) agree well with the simple calculation of ^30/2 « 4. 

From Figure[6]we can calculate a sky-averaged signal to noise 
ratios, (p), for a detectable source. For the example case above of 
using the ETs for the search and an optimally oriented source (p) = 
26.7. We know that the signal amplitude that led to this was A = 
500 ns, so assuming that we want a source to have (p) ^ 10 for 
detection we can scale the signal amplitude appropriately giving an 
effective sky averaged amplitude required for detection of 



( ^ I 



187 ns. 



(8) 



To set a scale (which will effect the absolute values of our event 
rates estimated below, but will not effect the relative ratios between 
them), we will say that a source at a distance of 1 Gpc produces an 
amplitude in the pulsar residuals of 100 ns. The observed amplitude 
is directly proportional to the source distance D, so for detection 
we require that 



The ratio of signal-to-noise ratios between an ET search (using 



D 



100 ns 



1 Gpc J \ (A) ns 



(9) 
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(a) 




8 10 12 14 16 



(b) 




(c) 



Figure 6. Hammer projection sky maps giving the signal-to-noise ratios for a source located at different sky locations for the IPTA observations with: | (a) | the 
signal appearing in the ETs of all pulsars for an optimally oriented source; (b) the signal appearing in the PTs for pairs (or multiples) of pulsars for an optimally 
oriented source and using the loudest pair for a given sky location; and, (c) the signal appearing in the PTs for pairs (or multiples) of pulsars for the worst 
source orientation and using the quietest pair for a given sky location. 
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or more generically 



D 



100 ns 



(10) 

lGpC/ V AllS J V(P)thresh/ 

where A is the observed signal amplitude that gives a sky-averaged 
signal-to-noise ratio of (p)a, and (p) thresh is the signal-to-noise 
ratio threshold for detection (so given the case in Equation [8] 
A — 500 ns, (p)a — 26.7 and (p)thresh = 10). This gives us a 
sky-averaged range, or horizon distance (similar to that used as a 
figure of merit for compact bin ary coalescences in ground-based 
gravitational wave detectors e.g. I Abbott et"aDl2005h , as a function 
of the observed amplitude. 

From Figure [6l we can calculate that for the ET search the best 
and worst case sky-averaged signal-to-noise ratios, (p) are 26.7 and 
6.4 respectively, giving horizon distances of 0.53 and 0.13 Gpc. For 
the PT search the best and worst case sky-averaged signal-to-noise 
ratios, (p) are 6.6 and 1.3 respectively, giving horizon distances of 
0.13 and 0.03 Gpc. 



3.2.2 Effective observation times 

For the ET search each pulsar observes the whole sky, so the ef- 
fective total observation time, T to t, is just the span of the residual 
observations T to t,ET — T res , but for the PT search the calcula- 
tion is more complex. For a pulsar pair the overlapping observation 
time that could contain a signal (assuming that the signal width is 
small compared to the data span) will be T res — At, where At is 
the maximum delay between PTs. Also, whereas in the ET search 
each pulsar observes the whole sky, in the PT search each pulsar 
pair will only have a fractional sky coverage that gives PTs within 
At (the ring on the sky in Figure [2]). So, for a particular pulsar pair 
the observation time factored by the sky coverage will be 



(T res - At) X fi, 



(11) 



where fi is the fractional sky coverage for that pair. The effective 
total observation time will therefore be T to t,PT = J^f T OD s,i, 
where N is the number of pulsar pairs. As we are assuming T res 
and At are the same for all pulsar pairs this becomes 



Horizon Distance = 0.53 Gpc 
Horizon Distance = 0.13 Gpc 
Horizon Distance = 0.03 Gpc 




(T res - At) X^/i. 



(12) 



Effective Observation Time (years) 



Figure 7. The source event rate as a function of effective observation times 
required to observed one event for a variety of detector horizon distances. 
The vertical dashed lines show the effective observation times for a search 
using the PTs (left) and ETs (right) with the IPTA pulsars. 



sources that would produce a 100 ns amplitude signal if at 1 Gpc) 
against total observation time if we want to be able to observe one 
signal. This is shown in Figure [7] along with lines corresponding 
to the effective observation times of an ET search and a PT search. 
We find that for the ET search, with an effective observation time 
of T tc 
and ^ 

assuming all sources are either in the best case or worst case ori- 
entations respectively. For a search using the PTs, with an effec- 
tive observation time of T to t = 11.5 years, source event rates of 
~ 9.6 Gpc -3 yr -1 and ~ 790 Gpc -3 yr -1 would be required to 
observe one event with the best case and worst case orientations 
respectively. This suggests that many sources should be observed 
in an ET search if you were to observed one in the PT search, or 
conversely if you observed one event in an ET search you would 
expect it to be very unlikely that a PT search would also produce 
an observed signal. 



= 20 years, source event rates of ~ 0.08 Gpc 3 yr 1 
5.5 Gpc -3 yr -1 would be required to observe one event 



So, for the PT search the lack of sky coverage for individual pulsar 
pairs can be compensated for by the effective increase in the tem- 
poral coverage i.e. each pulsar pair sees a different time epoch, but 
can only see that epoch for a small portion of the sky. 

For T res = 20 years and At = 10 years we find that for 
all IPTA pulsar pairs the sum of their fractional sky coverage 
is $^*=i fi — This means that T to t,PT = 20 years and 

T to t,ET = 11.5, which are comparable. In the next section we 
discuss what this indicates in terms of event rates. 



3.2.3 Source event rates 

If we assume that the rate of events that would produce an observed 
signal amplitude of 100 ns if at a distance of 1 Gpc, R (per unit 
volume per unit time), is the same throughout the Universe we can 
approximate the observed number of events, O, as 



Q = Rx 



4 o 

r D x 



(13) 



Using the horizon distances calculated above we can use 
Equation [13] to calculate the rate of these events (i.e. the rate of 



3.2.4 An expanded array 

Using the example of the IPTA above we can scale values to see 
how ET searches and PT searches would compare for a larger array 
of N pulsars given a fixed source event rate R. For an ET search, 
given that the observation time is fixed, the number of events that 
will be observed is dependent on the observing volume, which it- 
self depends in the horizon distance, which corresponds to the ob- 
served sky-averaged signal-to-noise ratio. The signal-to-noise ratio 
increases with the square root of the number of pulsars in the array, 
so overall the observed number of events becomes 



Oet 



4 



Deta 



1/2 



x T tc 



(14) 



where Det,ipta and Tet, tot are the horizon distance and effec- 
tive observation times for the ET search using the IPTA above. For 
the PT search the signal-to-noise ratio, and hence horizon distance 
and observed volume, is fixed (assuming that we still only use pairs 
of pulsars, rather than more), but the effective total observation time 
changes. If each pulsar pair sees on average a fraction / of the sky, 
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Number of Pulsars in Array 

Figure 8. The ratio of the number of events observable in an ET search 
compared to that observable in an PT search as a function of the number of 
pulsars in a PTA. 

and the number of pulsar pairs given N pulsars is N( N— 1)/2, then 
the effective observation time is proportional to fN(N — 1)/2, and 
the observed number of events becomes 

Opt = Rx ^7rD 3 PT , 1PTA x T PT fN(N - l)/2, (15) 

where L>pt,ipta is the horizon distance for the PT search using 
the IPTA above and Tpt — T res — At = 10 years as above. It 
can be seen that Equation Q31 scales as N 3 ^ 2 whereas Equation [HI 
scales as N 2 , so for a fixed rate (provided the sources and event 
rates are isotropic within the horizon distance) and a large enough 
array of pulsars the number of events observable in a PT search 
should overtake that observable in an ET search, i.e. the total ef- 
fective observation time available will compensate for the smaller 
range. 

Using the numbers for the IPTA (i.e. assuming a maximum 
delay between PTs of 10 years) we see that the average sky frac- 
tion observed by a pulsar pair is / = 1.15/435 = 0.0026. Using 
the horizon distances for the best case orientation for the ET and 
PT searches Figure [8] shows the ratio of the observable number of 
events in an ET search to a PT search as a function of the number 
of pulsars in an array. We see that to achieve equality in the number 
of events observed would require an array of ~ 400 000 pulsars. 
This is grea ter than the tot al expected number of active pulsars in 
the Galaxy (Lorimer 2008), so is not achievable. 



4 SUMMARY 

In searches for stochastic or burst sources of gravitational waves 
in PTA data the delayed PT is often ignored as it will generally 
be incoherent between residuals from different pulsars. We have 
looked at whether extra value could be gained when looking for 
short duration burst sources by looking for coherent signals in PTs 
from pairs, or more, of pulsars for specific source sky locations. 
In the era of the Low Frequency Array (LOFAR)ZSK A it could be 
possible to observe a few thousand millisecond pulsars JSmits et aD 
l2009h . with arou nd 100 of the best timed of these being usable as 
"arms" in a PTA ( Kramer 2012). However, our studies have mainly 
focused on the possibilities using the IPTA, containing 30 pulsars. 



We have shown that for a PTA timing array there can be sig- 
nificant areas of the sky containing a source for which the PTs of 
a signal can appear separated by relatively short times in pairs of 
pulsar residuals. In particular for the IPTA we find that this is the 
case for 67 per cent of the sky given maximum delays between PTs 
in two pulsars of 10 years. To accurately know which sky locations 
this is valid for would require precise (sub one per cent) knowledge 
of the PTA pulsar distances. This may be possible in the future with 
the SKA, but otherwise any search would have to take into account 
the distance uncertainties. For arrays of more pulsars the sky cov- 
erage becomes greater, e.g. for 90 pulsars almost the whole sky is 
covered, and there also become more source sky locations that give 
PTs observable in more than two pulsars. 

We have also compared the relative sky-averaged sensitivities 
of a search using the ETs, which are coherent between all pulsars 
in the PTA, with that of a search that uses the PTs in pairs of pul- 
sars. Considering an optimally oriented source, and using the pul- 
sar pair with the largest antenna pattern for a given sky location, we 
find that the ET search is on average about four times more sensi- 
tive. For the worst case source orientation, and using the pulsar pair 
with the smallest antenna pattern we find that the ET search is on 
average five times more sensitive. These are very comparable with 
what might be expected given that the ET search uses 30 pulsars, 
whereas the PT search is generally just using pairs of pulsar i.e. 
^30/2 « 4. An ET search would require a signal to occur within 
the observation window of the residuals, whereas each PT search 
would be looking back in time covering some time window depen- 
dent of the delay between pulsar terms. However, the ET search 
would cover the whole sky whereas the searches for PTs in pairs of 
pulsars would only cover small parts of the sky. For residual obser- 
vations covering 20 years and maximum delays between PTs of less 
than 10 years we have found that a PT search would have an effec- 
tive observation time of 11.5 years taking into account the limited 
sky coverage for each pulsar pair. We have converted these sensitiv- 
ities to horizon distances for a putative source, and along with the 
effective observation times, used them to estimate the event rates 
required to observe one event in the ETs and PTs of a simulated 
IPTA with 20 years of residuals. We have found that event rates for 
comparable sources show that you would have approximately two 
orders of magnitude higher chance of observing a signal in the ETs 
than the PTs for an all-sky search. However, this does not mean 
that PT searches are pointless. They could be complementary to ET 
searches in helping add more observation time to specific targets 
in the sky like large galaxy clusters. We have also assumed that all 
pulsar residuals have equivalent noise levels, but it could be that the 
least noisy pulsars could dominate the sensitivity and improve the 
comparisons. 

For arrays of larger numbers of pulsars we have seen that the 
scaling of the total effective observation time for a PT search rises 
faster than the search volume of the ET search. However, unfor- 
tunately the potential observable number of events in a PT search 
does not rise above that for an ET search until ~ 400 000 pulsars 
are in the array. Unfortunately such a large array is a completely 
unachievable goal. However, in larger arrays much more of the sky 
would be covered be cases when more that two pulsars have PTs 
separated by less than 10 years, so this would potentially bring that 
ratio down to a more hopeful level. 

In future work we plan to: investigate a more realistic timing 
array with more physical pulsar residuals; test and characterise a 
search routine; and, if possible, apply it to real data, or future IPTA 
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mock data challenge^. These studies will show whether the results 
here are more pessimistic than necessary. 
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APPENDIX A: THE ANTENNA PATTERN 

Given k is a unit vector pointing along the wave propagation di- 
rection from the source to the Earth we can construct the following 
basis vectors representing the wave propagation frame, 

k = ( — sin s cos (j)s 5 — sin s sin (f) s , — cos $ ) 5 (Al) 
I = (sin cf) s , cos (j)s i 0) , (A2) 
rri = (cos S cos (j) s , cos S sin cj) s , — sin S ) , (A3) 

and define the unit vector from the Earth (or more cor- 
rectly the solar system barycentre) to the pulsar n = 
(sin P cos P , sin P sin P , cos P ), where <j) SiP give the right as- 
censions and Q s ,p — tt/2 — declination for the source and pulsar 
respectively. From these the polarisation basis tensors can be de- 
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Figure Al. The absolute values of the 'plus' and 'cross' polarisation 
antenna pattern for pulsar J2129-5721 (a = 21 h 29 m 22f7, S = 
— 57° 21' 14" 1) over the whole sky. The pulsar location is marked by the 
black cross. 
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Figure A2. The absolute values of the 'plus' and 'cross' polarisation 
antenna pattern for pulsar J1603-7202 (a = 16 h 03 m 35!6, 8 = 
— 72°02 / 32' / 6) over the whole sky. The pulsar location is marked by the 
black cross. 



fined by 



ab id ib a b 

e i = I I — m m 



ab 



(A4) 

e x = l a m b +m a l b , (A5) 
and, with the inclusion of the polarisation angle ip, by 

e a + — cos (2'0)e+ b + sin (2 / 0)e' x b , (A6) 



- sin (2 , 0)e" b + cos (2ili)e% 



(A7) 



T he 'plus' and 'cross' ante nna patterns (e.g. see lAnholm et al.l2009l 
or lFinn & Lomm en 2010) are therefore given by 



n a n b 



1 ab 

z 1 + k • n 



(A8) 



Examples of these full sky antenna patterns for a couple of pulsars, 
with ip = 0, are given in Figures lATI and lA2l 

As we are not in the long wavelength approximation regime 
pulsar timing residuals will consist of components from the ET and 
PT, such that 



r(T) = {A + (t e )-A + {e))F + ^Xn) + 
(A x (t e )-A x (t p ))F x (^,kh) 



(A9) 



where, as we are dealing with the PTs in this paper, we set the time 
of the PT as t p = T — U, where T is the pulsar proper time at 
the solar system barycentre and ti is time at the observed midpoint 
of the burst in the pulsar term (as seen for a particular pulsar), and 
therefore the CTtime is t e = T — U + (1 + k • hi)di/c. 

When creating a coherent signal in two pulsars' residuals 
the signal model for the second pulsar will have to be shifted by 
Atij = (1 + k • hi)di/c — (1 + k • hj)dj/c with respect to the 



first pulsar i.e. if we have n (T) then to coherently combine it with 
another dataset that would have to have r2 (T+ At 12), with U = t\ 
for both models. 



